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The dynamics of expansion of the Universe and evolution of scalar perturbations are discussed 
for the quintessential scalar fields Q with the classical Lagrangian L = ^Q-iQ' 1 — U(Q) satisfying 
the additional condition w = const or c\ = 0. Both quintessential fields are studied for the same 
cosmological model. It is shown that the accelerated expansion of the Universe is caused by the 
effect of rolling down of the field to minimum. At the early epoch the contribution to dynamics of 
the quintessence with w = const is negligible (like that of cosmological constant) while quintessence 
with c\ — mimics dust matter. In future the scalar field with c\ = will mimic cosmological 
constant. 

The systems of evolution equations for gauge- invariant perturbations of metric, matter and 
quintessence have been analysed analyticaly for the early stage of the Universe life and numerically 
up to the present epoch. It is shown that amplitudes of the adiabatic matter density perturbations 
grow similarly in both models (and like in ACDM-model) , but time dependences of different ampli- 
tudes of the quintessence perturbations are varied: gauge-invariant variables and Ds decay 
from initial constant value after the particle horizon entry while and grow at the early 
stage before the horizon entry and decay after that - in the quintessence-dominated epoch, when 
gravitational potential starts to decay - so, that at the current epoch they are approximately two 
orders lower than matter ones on the supercluster scales. Therefore, on the subhorizon scales the 
quintessential scalar fields are smoothed out while the matter clusters. 

It is also shown that both quintessential scalar fields suppress the growth of matter density 
perturbations and the amplitude of gravitational potential. In these QCDM-models - unlike ACDM 
ones - such suppression is scale dependent and more visible for the quintessence with c\ = 0. 
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I. INTRODUCTION 



Cosmological observations of the last decade surely 
assert that the main part of the energy density of the 
Universe - more than 70% - belongs to the unknown 
essence, called "dark energy". Its cosmological mission 
is to provide the accelerated expansion of the Universe, 
revealed from exploration of SN la's in the distant galax- 
ies and temperature fluctuation power spectrum of cos- 
mic microwave background. The cosmological ACDM- 
model, based on the Einstein equations with cosmologi- 
cal constant (see 0, [H, |43[ and references therein) , de- 
scribes very well almost whole set of the observational 
data on dynamics of expansion of the Universe and for- 
mation of its large-scale structure. But physical interpre- 
tation of the cosmological constant is rather problematic 
0, 0, H HE US El • Therefore alternative approaches - 
new physical fields (classical scalar field - quintessence, 
tachyon field, k-essence, phantom field, quintom field), 
Chaplygin gas, gravity and general relativity modifica- 
tions, multidimensional gravity, branes and others - are 

intensively analysed (see reviews Silo, mm mug 
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and special issue of Gen. Relativ. Gravit., 2008, v. 40) 
now. Up to now none of them has crucial preferability 
from observational or theoretical point of view. Therefore 
each of them must be comprehensively studied. Here we 
restrict ourselves to quintessential scalar fields with clas- 
sical Lagrangian L — Q-iQ'^/2 — U(Q) in the dark energy 
- matter dominated Universe. 

The quintessence model can be defined by setting of 
the appropriate potential U (Q) or equation of state (EoS) 
parameter wq = Pq/c 2 pq. There is a dozen or more 
physically-motivated shapes of the potential U(Q): ex- 
ponential, double exponential, exponential with inverse 
power, power-law, etc. The dynamics of such scalar fields 
is intensively studied (see review [ll|). The EoS param- 
eter of dark energy completely defines the background 
dynamics as well as the evolution of cosmological per- 
turbations [U [12, Hl| . Since observational data on SN 
la magnitude - redshift relation and cosmic microwave 
anisotropy give relatively narrow ranges of dark energy 
density and EoS parameter values, it looks quite attrac- 
tive to establish the potential U(Q) using these data and 
analyse the background dynamics and perturbative prop- 
erties of such scalar field which are not studied widely 
enough. 

In our previous papers we have constructed the po- 
tentials of scalar fields with classical and tachyonic La- 
grangian leading to the constant EoS parameter wq — 
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const [40( and analysed the background dynamics and 
perturbative properties of such scalar fields [4l[. It was 
shown that cosmological model with cold dark matter 
and such types of the scalar field (QCD A/- model) agrees 
slightly better with the accessible today observable data 
than the ACDM-model. But difference of quantitative 
merits of goodness is not large enough to pick out one of 
them at confidential level of la. Since the degeneracies 
between model parameters of dark energy and cosmolog- 
ical parameters [H H3, HE HI S3 exist for the back- 
ground dynamics, the complete analysis of linear den- 
sity perturbations in both dark matter and dark energy 
components is important for improvement of dark en- 
ergy observational tests. Among the large number of free 
quintessence parameters and unknown initial values of 
quintessence perturbation modes there is only small part 
of models, for which the evolution of perturbations has 
been studied. The general conclusion is that magnitudes 
of dark energy density perturbations on scales smaller 
than horizon are essentially lower than corresponding 
magnitudes of matter density ones. But character of their 
evolution depends strongly on the scalar field model (its 
potential, time variation of EoS parameter, sound speed, 
etc.), initial conditions, scale of perturbations and gauge 
(see for example |, |, [H, 0, Q E3, US H El ) ■ 



Here the special attention should be paid to the EoS 
parameter of dark energy wq, which can be constant or 
varying in time. The temporal variation of the dark en- 
ergy EoS parameter is often presented by linear fitting 
formula with two or three [24| parameters to be es- 
timated. Other functional dependences of wp on scale 
factor or redshift can be found in [TJ [27], Hf|. Here 
we study the parametrization of the equation of state, 
which needs only 1 additional quantity with clear physi- 
cal meaning - the adiabatic speed of sound c 2 = pq j c 2 pq 
(the analysis of generalized dark sector components can 
be found in the early works [5l],[23]). ^ n general, c 2 is the 
unknown function of time. However, taking into account 
the simplicity we restrict ourselves to c 2 = const, so it is 
regarded only as the second physical parameter defining 
the equation of state of dark energy (the first one - the 
present value of wq). 



In this paper we undertake the comparative analysis 
the evolution of gauge-invariant variables of the scalar 
perturbations in the model with non-relativistic matter 
(pm c 2 pm) and scalar field which we define by classi- 
cal Lagrangian with potential constructed for two cases 
(wq — const and c 2 = 0) in the concordance cosmologi- 
cal models. These cases have been chosen because they 
allow us to obtain analytical solutions, which seems to 
look very attractive in the world of numerical compu- 
tations. We assume the adiabatic initial conditions for 
matter and dark energy scalar perturbations. 



II. BACKGROUND COSMOLOGICAL AND 
SCALAR FIELD MODELS 

We consider the homogeneous and isotropic flat Uni- 
verse with metric of 4-space 

ds 2 = g lJ dx l dx j = c 2 dt 2 - a 2 (t)5 a 0dx a dx 13 
= a 2 (n)(dn 2 - S aP dx a dx 13 ), 

where the factor a(t) is the scale factor, normalized to 
1 at the current epoch to, 77 is conformal time (cdt = 
a(ji)dn). Henceforth we also put c = 1, so the time vari- 
able t = xq has the dimension of a length. Here and 
below the latin indices i, j, ... run from to 3, the greek 
ones - over the spatial part of the metric: v, p, ...=1, 2, 
3. 

If the Universe is filled with non-relativistic matter 
(cold dark matter and baryons) and quintessence which 
interact only gravitationally (minimal coupling) then the 
dynamics of its expansion is completely described by the 
Einstein equations 

%-^ = 8 7 rG(i;f ) +^ Q) ), (1) 

where Rij is the Ricci tensor and , 2y - energy- 
momentum tensors of Matter (M) and Quintessence (Q). 
If these components interact only gravitationally then 
each of them satisfy the differential energy-momentum 
conservation law separately: 

T i(M,Q) = Q (2) 

(here and below ";" denotes the covariant derivative with 
respect to the coordinate with given index in the space 
with metric gtj). For the perfect fluid with density P(m,q) 
and pressure P(m,q)> related by the equation of state 
P(M,Q) = W( M ,Q)P(M,Q), it gives 

P(M,Q) = ~3-P(Af,Q) (1 + W (M ,Q) ) (3) 

(here and below a dot over the variable denotes the 
derivative with respect to the conformal time: "'" = 
d/dn). The matter is considered to be non-relativistic, 

so w M = and p M = p$a~ 3 (here and below "0" de- 
notes the present values). 

We assume the quintessence to be a scalar field Q(x, n) 
with classical Lagrangian 

L= l -Q^ l -U{Q), (4) 

where U(Q) is the field potential. We suppose also the 
background scalar field to be homogeneous (Q(x, 77) = 
Q{n)), so its energy density and pressure depend only on 
time: 

PQ(V) = ^Q 2 + U(Q), PQ (n) = ^Q 2 - U(Q). (5) 
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Then the conservation law ([2]) gives the scalar field evo- 
lution equation (called the Klein-Gordon one) 

Q + 2aHQ + a 2 — = 0, 
dQ 

where H — a/ a 2 is the Hubble parameter for any moment 
of conformal time n. 

We specify the model of quintessence using two ther- 
modynamical parameters: the EoS parameter wq = 
Pq/pq an d the adiabatic speed of sound c 2 = Pq/pq- 
In general case they are connected by equation 

dw /d In a 2 

M =w " Ca 

(here and below we omit index Q for wq). If the time 
dependence of uu is known then c 2 is defined unambigu- 
ously, if c 2 is determined then the initial value Wq must 
be defined additionaly, so, the EoS parameter has 2 de- 
grees of freedom: a function and a constant. For other 
parametrizations see [III H3, [H[ . Since the constraints 
for time dependence of w or c 2 are not established well we 
consider two simple cases: w = const and c 2 = const. In 
the first case c 2 = w and in the second one l+w(a) = (1 + 

c 2 a )(l+w )/ (l + w - (w Q - c 2 a )a 3< - 1+c ^ . This equation 

has obvious asymptotical behaviour: when a — > w — ► c a 
and when a — > 00 w — ► — 1. So, at early epoch the dark 



energy mimics dust matter (w w 0) for c 2 — or radia- 
tion (w « 1/3) for c 2 = 1/3. In future such scalar field 
will mimic cosmological constant (w ~ —1). The time 
dependences of EoS parameter for both cases are shown 
in FigQ] The equation ([3]) has the analytical solutions 
for two cases: 



w = const: pq(a) = Pq 



-3(l+w) and 



• c; 



= 0: p Q {a) = p ( Q ] [(1 + Wo)a~ 



so it's possible to simplify formulae and calculations and 
we will analyse only this two cases now. 

If the parametrization of EoS parameter is given, it is 
possible to apply reverse engineering and construct the 
fields Q and potentials U(Q). From (JSJ) one simply ob- 
tains: 



Q(o) - Qo = ± 



a Vpq(1 + w) 
aH 



U(a) 



pqO- 



If the integral for Q can be expressed via functions that 
could be inverted to obtain a(Q — Qo), then U(Q — Qo) 
can be easily written in analytical form. 

So, from the Einstein and field equations we deduce 
the time dependences of the Hubble H and acceleration 
q parameters as well as the evolution of the scalar field 
Q and potential U(Q): 



H = Hoa 
Q(a) - Qc 



± 



1 VI + w 



ll-fi c 



q= 2 



(1 + 3w)Qqa 



-3w 



l-Slc 



-Sw 



In 



V(i - n Q )a 3 ™ + n Q - yn^ 1 + y/Uc 



y (1 - n Q )a 3w + n Q + 1 

1 



ch V6^G(Q - Qo) 



sh { V6nG(Q - Qo) 



1 2- 



(6) 
(7) 

(8) 



for w = const and 



H = H a 2 \/l + Qqw - QqWoo, 3 , 



11 + WoQq + 2wo^Iqcl 3 
2 1 + Oqw - flQW a 3 



Q(a) - Qo = ± 



U(Q-Qo) = 



2V6^G V 

3H 2 n Q (i 



8vrG 



n Q (l + w ) I yfi + Vi Q wo{l - a 3 ) - yT+ » Q w 1 + y 7 ! + Oqwp \ 
1 + Q Q w n ^ ^1 + ^(5^0(1 - a 3 ) + a/1 + f^QU'o 1 - V 1 + ^QWo / 

cfo ( v/^rGCQ-Qo). 



(9) 
(10) 



w ) 



' 1 + ^qwq 

Q Q (1 + w ) 



4 



3/1 V6^G(Q-Q ) 1 



' 1 + OqWq 

Qq(1 + wo) 



-, 2 



8ttG 



(11) 



for c 2 = 0. For both models there are 2 independent 
solutions for the field (the growing one corresponds to 
sign "+" and the decaying one to sign "-") and 2 sym- 
metrical with respect to Q — Qq potentials exist. How- 
ever, the physical consequences of both these solutions 
are the same (40l | , so from now we restrict ourselves only 
to the growing one. The variety of scalar field poten- 
tials was presented in 11111 . The potential for w — const 
can be also found in [3a, 53], the potential for c 2 = 
belongs generally to the family of double exponential po- 
tentials (with additional constant term), but both they 
differ from the physically-motivated ones, for which the 
evolution of scalar linear perturbations was studied by 
other authors. 

We must note that the asymptotic behaviour at a — > 
of the expansion rate H(a) and acceleration parame- 
ter q(a) in both cases is the same and similar to that 
in ACDM- model: H cx a" 3 / 2 , q -> 1/2. At current 
epoch the parameters of expansion dynamics are the 
same (Ho and qo = (1 + 3wQ,q)/2) for both models. 
But their asymptotic behaviour at a — > oo is differ- 
in w — const quintessence H — » Hoa~^^ 1+w ^ ^Qq, 
(1 + 3w)/2, pq — ► and in c 2 — quintessence 

7T~ 1 (°) 

WqHq, q -> -1, Pq -> -W p Q - 



cnt 

q - 

H 



3w /(l — Qq) always while in the c 2 = 
Pq/pm -> (1 + io )O q /(1 - n Q ) 



ffo^uiofig, g -> -1, pq -WqPq 1 . The energy 
densities of both fields evolve similarly but have different 
asymptotic regimes: in the quintessence with w = const 
Pq/pm = Qqci 
quintessence at a 

and at a — > oo Pq/pm — * — wo^Qa 3 /(l — fig). So, the 
scalar field with c 2 = behaves as cold dark matter at 
the early epoch and will mimic the cosmological constant 
in far future. 

The different asymptotic behaviour of these fields is 
caused by their intrinsic properties. In the w = const 
quintessence the negative pressure stiffly follows its en- 
ergy density and their relation is always constant. In 
the c 2 = quintessence the negative pressure is always 
constant: pq — 3HqWoQq/8ttG. So, it is insignificant 
in the early epoch when a — ► and pq — ► oo for the 
model of the Universe filled only with dust matter and 
quintessential dark energy, and important in the late one 
when w — > — 1. 

The dynamics of expansion of homogeneous Universe 
in the model with non-relativistic matter and quintessen- 
tial scalar field with c 2 = and w — c 2 — const is 
shown in Fig[TJ For both models we assume best fit- 
ting cosmological parameters from [43| (Qq = 0.745, 
w = w = -0.915, fl M = 0.255, h = 0.7). For com- 
parison we show also the corresponding dependences in 
ACDM-model with fl A = 0.74, Q M = 0.26 and h = 0.73 



am. 

We have constructed the pot entials of quintessential 
scalar fields with w — const [4(| and c 2 = for QCDM 
cosmological model with best fitting parameters obtained 
from WMAP and SNIa data H|. The evolution of fields 
Q(a), potentials U(a) and rolling down of the fields Q to 
the minimum which is located at Q — > oo (a — > oo) are 
shown in FigfSJ The discussion of influence of parameter 
determination uncertainties on potential of field with w = 
const can be found in [401 ] . 

So, the difference of the homogeneous Universe expan- 
sion dynamics in ACDM- and such QCDM-models is too 
small to discriminate them using the avialable datasetets. 
That's why in the next sections we will analyse the lin- 
ear stage of growth of scalar perturbations of matter and 
dark energy. For this we will use gauge-invariant ap- 
proach developed by 0, 



III. EVOLUTION OF SCALAR 
PERTURBATIONS 



For analysis of the scalar linear perturbations the 
conformal-Newtonian gauge with space-time metric 



ds 2 = a 2 (ri)[(l + 2*(x,7 ? ))d77 2 
-(1 + 2^(x,n))S aP dx a dx 13 ] 



(12) 



is convenient. Here ^(x, 77) and $(x, 77) are gauge- 
invariant perturbations of metric |2| called Bardeen's po- 
tentials. If proper anisotropy of medium equals zero 
then \&(x, 77) = — $(x, 77). Dust matter and scalar 
fields have such property [23j. In the linear per- 
turbation theory the Fourier decomposition is used, 
so spatial dependences of all variables can be substi- 
tuded by corresponding Fourier amplitudes. For ex- 
ample, ^(x, 77) — > ^(k, rj), where k is wave number. 
Henceforth, saying about metric ^(x, rj), matter density 
<5 (M) (x, 77) = (pm(x, 77) - Pm(v))/pm(v), its peculiar ve- 
locity y( M )(x,7y), scalar field SQ(x,j]) = Q(x,rj) - Q(rj), 
its energy density perturbations 5^ (x, rj) = (pg(x, 77) — 
PQiv))/ Pq( 7 i) etc we mean their Fourier amplitudes 
S {M) (k,r]), (*,»?), 8Q(k,rj), 5^(k,rj), etc. 
The metric (^(k, rj)), matter density and velocity per- 
turbations (5( M '(k, 77), V( M ^ (k,n)) as well as scalar field 
perturbations (5Q(k,rj), S^'(k,rj), V^'(k,rj)) in the 
conformal-Newtonian gauge are gauge- invariant variables 
[23| . The energy density and velocity perturbations of 
quintessence, 5^ and 7^' , are connected with the per- 
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FIG. 1: Top: the dependence of EoS parameter w on scale fac- 
tor a for c\ = 1/3, 0, —1/3 and w — c\ = const. Middle: the 
dynamics of expansion of the homogeneous Universe in the 
model with non-relativistic matter and quintessential scalar 
field with c 2 a = and w = c 2 a = const with best fitting cosmo- 
logical parameters from Spergel et al. (2007) (fig = 0.745, 
w — —0.915, Qm = 0.255, h — 0.7): matter and quintessence 
densities in units of the critical one. Bottom: the evolution of 
acceleration parameter. For comparison we show also the cor- 
responding dependences for ACDM-model with Qq = 0.74, 
Q.M = 0.26 and h = 0.73 (Spergel et al. (2007), Apunevych 
et al. (2007)). 



turbation of field variable SQ in following way: 

S^> = (1 + W - * H r-^ — 

\Q Q 2 dQ 



V«3) = 



kSQ 



Other non-vanishing gauge-invariant perturbations of 
scalar field are isotropic pressure perturbation 



(q) _ 1 + w / 5Q_ 

Xh «> V Q 



a 2 SQ dU 



and intrinsic entropy 
r (Q) 



J<2) 



^6 {Q) . 
w 



The density perturbation of any component in the 
conformal-Newtonian gauge D s = 5, which is gauge- 
invariant variable, is related to the other gauge-invariant 
variables of density perturbations D and D g as: 



D 



D g +3(l+w) 



aV 
a k 



D a +3(l+w)~, (13) 
a k 



where D s , D, D g and V correspond to either M- or 
Q-component. Here D g is the density perturbation in 
the rest frame in which the fluctuations of the curvature 
scalar of the constant time hypersurface vanish and D 
corresponds to the rest frame in which the 4-velocity is 
orthogonal to constant time hypersurface [23| . 

The intrinsic entropy of quintessence can be pre- 
sented via gauge-invariant Q-perturbations as follows: 



wT 



(Q) 



a-ci)D 



(Q) 



(14) 



This equation shows that the intrinsic entropy for scalar 
perturbations of quintessence with c 2 ^ 1 is non-zero 
when proper energy density perturbation (mea- 
sured in synchronous comoving gauge) of quintessence 
is non-vanishing. In the first case (w — const) wT*®' — 
(l-w)D (Q \ in the second one (c\ = 0) wT^ = D®h In 
the case of perturbed quintessence dissipative processes 
generate entropic perturbations, so we have the sound 
speed c 2 denned by more general relation: c 2 = Spq/Spq. 
The intrinsic entropy perturbation can be presented in 
the form: wT^ = {c 2 s - c 2 a )D^ [23|L_For the scalar 
fields with classical Lagrangian c 2 = 1 [18|, [45[ . 



A. Evolution equations 

Evolution equation for scalar field perturbation 
5Q(k,rj) can be obtained either from Lagrange-Euler 
equation or from energy-momentum conservation law 



T, 



(Q) 



0;/ 



0: 

SQ + 2aH5Q 

-4*g = o. 



,d 2 U 
dQ 2 



S Q + 2a 2 ^ 



dQ 



(15) 



Thus, evolution of quintessence perturbation depends on 
field model (U(Q)), gravitational potential 'J, expansion 
rate of the Universe H and scale of perturbation k. 

The linearised Einstein equations for gauge-invariant 
perturbations of metric and energy-momentum tensor 
components are 
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FIG. 2: Top: the evolution of fields Q(a) (left), potentials U(a) (right) for quintessence with w = const and c\ = 0. Bottom: 
rolling down of the fields Q to the minimum U(Q) = which is located at Q — > oo (a — + oo) for cases of w = const (left) and 
c 2 a = (right). 



£>'«> + *(l -«,)£>(«)- 
y a * 

y,{Q) _ ly{Q) _ 

a a 2 H 

T , * 4nG ( 

* + TT [PM— 7— 

a H \ k 



a 2 H 
k D 



9 (1 f a)H \ (1 + w)V™ + 9(1 + w)(l - cl)- = 0, 
K J a 



(Q) 



o?H 1 + w 



= 0, 



p Q (l + w) 



0. 



(16) 

(17) 
(18) 



Here and below a prime denotes the derivative with re- 
spect to the scale factor a. The conservation equations 
for matter density and velocity perturbations 5Tj.^ M ^ — 

in terms of the gauge-invariant variables D g M ' and 
y(M) arc f n owm g : 



(19) 
(20) 



They are connected with the dark energy ones only via 
^ and are the same for both models of quintessence. 



So, in each case we have the system of 5 first- 
order ordinary differential equations for 5 unknown func- 
tions *(fe,o), D ( g M \k,a), V<> M \k,a), D ( g Q) (k,a) and 
V(Q\k,a). From this systems of equations it is easy to 
obtain the systems of 2 second-order ordinary differential 
equations for 2 unknown functions ^(k,a) and SQ(k,a): 



-.-tc^f N / facna( i + „) ii °*y + ^-"o . o, (2D 
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SQ" + 



_ -|(l-H») 



An 



H V 8ttG 



a 



9(1 



w)(2 + w) 9w(l 



3(1 - w)* 



4a 2 



= 



4a 2 



i£) nga -3(l-H»)|| 



<5g 



(22) 



for ui = const or 



2a6Q' + 3SQ 



2a 



= 0, (23) 



5 3 n H$\6Q' ( k 2 9 9 H 2 \ ^ _sH Q I 3 " " " 



4a*' + 3* 



(24) 



for eg = 0. 

Using their solutions (four fundamental) for *(/c,a) 
and Q{k,a) and the constraint equation 

- = iirGa 2 (p M D {M) + PqDM) , (25) 

it is possible to find the values of D g M ^ (k, a), V^ M ^ (k, a), 

D g Q \k,a), V( Q \k,a) and T^)(k,a). The equations 
(flT? ]) -(|2T) ]l can be substituted by one second-order equa- 
tion 

i^W + ( 2- g )_2_ + * * = . (26) 
J a aril* 

The systems of equations ([2"T ]) -([2"2" j) and (|2"3 > p3 1) de- 
scribe the evolution of perturbations of gravitational * 
and quintesscntional SQ fields and their coupling. A few 
important conclusions can be deduced from qualitative 
analysis of these systems: 

• The coupling of *— and SQ— field is modulated 
by the value 1k /r2q(l + w) for w — const and 
^/£!q(1 + w ) for c 2 =0. So, if w — w = — 1 
then both fields evolve independently. Since obser- 
vational data prefer current w close to —1, so their 
coupling is weak. 

• Evolution of SQ— field depends on relation of scale 
of perturbation to horizon explicitly while the de- 
pendence of *— field is implicit (through the latter). 

• The system of equations for u;-quintessence (|21[1 - 
(|2"2")l allows the asymptotic behaviour * — ► const, 



SQ — > const when a — > oo with relation between 
them SQ = */^6vrG(l + w). 

• The system of equations for c 2 -quintessence 
()23 jl - ([24|) allows the asymptotic behaviour 

* — > const, SQ — ► const when a — ► 
with initial relation between them <5Q = 
v/no(l + u; )/67rG(l + figwo)*- 

• From equation (f2"6"|) it follows that Dg M ' 1 » const 
for superhorizon perturbations (fc <C a 2 H). If 

* = const and g = 1/2 at a ^ then either 
Z)( M) = const -2ak 2 y/3H${l-n Q ) for to = const 

or Z^ M) = const - 2afc 2 */3# ( 2 (l + n Q w ) for 
c 2 = and for * < it begins to grow slowly from 
the constant value. The decay of 4" and transition 
from deceleration to acceleration slow the growth 

f4 M ». 

B. Initial conditions 

Analysis of the background dynamics presented in the 
previous section has shown that both QCDM-models are 
matter-dominated in the early Universe (Fig. [I}. In the 
QCDM-model with w — const the ratio pm/pq — ► oo 
when a — > 0, while in the QCDM-model with c 2 = 
Pm/pq — ► (1 — ^q)/(1 + wo)r2Q and iu — * when a — > 0. 
The adiabatic growing mode of perturbation in the non- 
relativistic matter-dominated Universe can be specified 
by the condition * = const (* = 0). Adiabaticity 
condition in two-component model (Sm-.q = D g M ^ — 
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Df ] /(l+w) = [3) gives D ( g M) = D { p ] /{l+w). These 
conditions, constraint equations written for hypersurfacc 
Vinit "C rj Q (a init <1) and the analytic asymptotic solu- 
tions (see next subsection) lead to the following adiabatic 
initial conditions: 



{ ry\ 2 k ^ init 

v init — 77 T7 t ~ r< == V ^in-it ; 



3 i?o a/1 — fig 
-5(l + tu)*imi, 



D (Q) = 

9 init 

V (M) _ 2 k ^init r- 

» in.it ~ ^ -r V u>imt : 



3 H — fig 



9 init 

for u; = const and 



^^^^ ?"rv,?'f — 



2 fc 



D (0) 

9 init 
^ init 



D {M) 
9 init 



3 H o yj\ + SIqWq 
-5^ init, 

2 fc ^init 



3 o a/1 + Qqw 

— d^init 



y/Oi- 



\J Q*ini 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 
(33) 
(34) 



for c\ = 0. 

Therefore, the growing mode of adiabatic pertur- 
bations in two-component (non-relativistic matter and 
quintessence) medium is defined by single value - initial 
gravitational potential ^i n it- 

Since the non-adiabatic initial perturbations are 



strongly constrained by WMAP data, in this paper we 
restrict ourselves only to adiabatic initial conditions. 



C. Asymptotic and numerical solutions 

In order to analyse the evolution of gauge-invariant 
variables of matter and quintessence perturbations we 
must solve the system of equations (fl"6|) - (fT8)) together 
with (Hnj)-(Eni) numerically for initial conditions (j2"7]) - (|3T))) 
or (f3"Tj) - (|3"4")) respectively. But before we propose the anal- 
ysis of these systems of equations in the early epoch 
(a <C 1), for which the analytical solutions are known. 
So, the system of equations (f23|)-([24l for a <C 1 can be 
simplified as 



„ 7*' 3n Q (l + w )-f 6TrGn Q (l + w ) 



2 a 2 1 + QqWq 



1 + Qqwo 



2aSQ' + 35Q 
X 2^ _0 ' 



3 n Q (l + w ) 



An „ 9_ 

V + 2 a + 2a 2 ° V y 8ttG 1 + Q Q w 

4a*' + 3* 
x - = 0. 



(35) 



(36) 



This system of equations has 4 fundamental solutions, so 
it is possible to write the general solution in the form: 



4< = 



SQ = 



02 



n Qmo +8n Q -7 

1 + "Q™0 



1+177 



1 



in Q (i + w ) 



\/67rGU l + fl Q w 



2 at fl Q (l + w ) 



1 + £Iqw j n Q w + 8S1q - 7 
1 + flqwo 



C 3 a 



i+n Q ™ 



CaCI 



The first two solutions, noted by the constants of inte- 
gration C\ and C*2, are well known growing and decay- 
ing modes of adiabatic perturbations in the dust matter- 
dominated Universe. The next two solutions, noted by 
the constants of integration C3 and C4, are due to possi- 
ble entropy initial conditions and intrinsic non- vanishing 
entropy of quintessence. Really, the condition = 
leads to 1 second-order equation which has two dust-like 
fundamental solutions: 



For the quintessence with w = const solutions for \t are 
the same and D { S Q) = -2(1 + w) [px - fCafl-iV 

The quantities Z^ M) (/c,a), V (M \k,a) and V^(k,a) 
can be found using the equations (fT3|) - ([2"0|) and (|25| . The 
relations between them are presented in the previous sub- 
section as the set of initial data (f27|) - (f34|) . 

We have integrated numerically the systems of equa- 
tions <pTH|) - (1 1 8 [) for w = const and for c\ = together 
with (TXni)-(Enj) for adiabatic initial conditions (12"7)) - (13"0|) 
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FIG. 3: The evolution of gauge-invariant amplitudes of perturbations in matter (top) and quintessence (bottom) for two 
models of quintessence: w = const (left column) and c\ — (right column). The corresponding scale of perturbations is 
k = O.OOlMpc" 1 and the cosmological parameters are Qq = 0.745, w = —0.915, Qm — 0.255, h = 0.7. 



and l[31"1 ) -(|34| ) using the publicly available code DVERK 1 . 
We assumed ai n u — 10~ 10 and integrated up to a = 1. 
The evolution of perturbations is scale dependent, so we 
performed calculations for k — 0.0001, 0.001, 0.01 and 
0.1 Mpc^ 1 for the cosmological model with parameters 
n Q = 0.745, w = -0.915, Q M = 0.255, h = 0.7. The 
evolution of gauge-invariant variables of matter perturba- 
tions D g M \ D { S M \ D( M \ for two scales k = 0.001 
and 0.01 Mpc -1 is shown in top panels of Fig[3]and FigfJ] 
In the bottom panels the analogical gauge-invariant vari- 
ables of quintessence perturbations (Dg, Ds , 

wT^Q}) are presented (for c\ = the curves wT^ 
and I)W overlap). The evolution of gauge- invariant 
gravitational potential \P is shown in all panels for com- 
parison. All plots are shown for the following range: 
0.001 < a < 1. 



IV. DISCUSSION 

From the top panels of FigO [4] it follows that the 
magnitudes of the adiabatic matter density perturbations 
grow similarly in both models (and like in ACDM-model), 
but time dependences of magnitudes of the adiabatic 
quintessence energy density perturbations are more var- 



1 It was created by T.E. Hull, W .H.Enright , K.R. Jackson in 19 76 
and is available at |http: //www.cs.toronto.edu/NA/dverk.f.gz 



ied (bottom panels of the same figures): gauge- invariant 
variables D g ^ and decay from initial constant value 
after the particle horizon entry while and 7'^ grow 
at early stage before the horizon entry and decay af- 
ter that - in the quintessence-dominated epoch, when 
gravitational potential starts to decay. The perturbation 
shown in Fig[3] enters the particle horizon (77(a) = ir/k) 
at a « 0.03 for w-quintessence and at a w 0.04 for c\- 
quintessence. The perturbation shown in Fig[5] enters 
the particle horizon at a w 0.0004 and a w 0.0005 for w- 
and c^-quintessence respectively. The particle horizon at 
current epoch (770) in the cosmological model with pa- 
rameters n Q = 0.745, w = -0.915, fl M = 0.255, h = 0.7 
and w-quintessence equals « 14970 Mpc. In the model 
with c^-quintessence it is « 13810 Mpc. At early epoch 
D {Q) 

oc a for both models of quintessence. After the par- 
ticle horizon entry the amplitudes start to decay slowly 
in the matter-dominated epoch and decay fast in the 
quintessence-dominated one. At asymptotic regime for 
quintessence model with c\ — approximately oc 
a~ 3 . In the quintessence model with w — const the 
transition epoch is extended in time. In FigO we show 
the dependences of ratios of quintessence density pertur- 
bations to matter density ones in conformal-Newtonian 

gauge (D^ I Dg^ ) on scale factor for perturbations with 
the scales k = 0.0001, 0.001, 0.01 and 0.1 Mpc' 1 . These 
curves emphasise the difference of evolution of perturba- 
tions in ordinary matter and quintessence as well as the 
similarity of behaviour of perturbations in two models 
of quintessence. The magnitudes of quintessence density 
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FIG. 4: The same as in FigEJfor scale k = O.OlMpc" 1 . 



perturbations in units of matter ones in both models at 
current epoch are close although their initial magnitudes 
differ by order. The magnitudes of quintessence density 
perturbations with scale less than particle horizon are 
lower than corresponding magnitudes of matter density 
perturbations by factor « (23000fc) 2 so, that for scale 
k = 0.01 D { S Q) /D { S M) « 2 x 10" 5 . Therefore, on subhori- 
zon scales the quintessential scalar fields are practically 
smoothed out while the matter clusters. 

The tests for choice of the type of dark energy are 
based on the results of its action on luminous matter 
and cosmic microwave background. So, the key question 
is how these types of quintessence affect the growth of 
matter density perturbations and the time variation of 
gravitational potential. From top panels of FigJ3] and [H 
we can see that they are more suppressed for c\ — than 
for w = const and for perturbations with smaller scale. 
In order to illustrate this effect in FigJS] we present ra- 
tios D^a in u/D^2a and */* init for scales k = 0.0001, 
0.001, 0.01 and 0.1 Mpc -1 . We can see that scale depen- 
dence of suppression of magnitude of matter density per- 
turbations as well as of gravitational potential is strong 
for c\ = quintessence and weak for w = const one. 
In the ACDM-model it is scale-independent @. (In the 
Einstein - de Sitter model both ratios are equal to 1 for 
all times and scales). These ratios are substantial for 
calculations of magnitude of the matter density power 
spectrum at different redshifts and the angular power 
spectrum of CMB temperature fluctuations in the range 
of scales of the late integrated Sachs- Wolfe effect. 

Evolution of quintessence perturbations depends on 
scalar field model (i.e. its Lagrangian and potential), 



contents of the Universe, coupling of the quintessence to 
other components, initial conditions and scale of pertur- 
bations @, i, GJ, [jj, [M [H, 13 . Here we have analysed 
the evolution of scalar matter and quintessence pertur- 
bations for potentials of scalar fields with classical La- 
grangian constructed to give either w = const or = 0. 
Therefore, obtained here results could be compared to 
the results of other authors only qualitatively. Evolution 
of EoS parameter in our c\ = 0-model (FigfT]) is simi- 
lar to that of [HI • Despite the other cosmological model 
and potential of scalar field, the qualitative behaviour 
of quintessence perturbations is close to obtained here: 
Dg is const when the perturbation is outside the par- 
ticle horizon and decays when it enters the horizon. The 
growth of magnitude of quintessence density perturba- 
tions long before the horizon entry in synchronous gauge 
was shown by [l2T ] (models with w — const in Fig. 3). 
Our results for evolution of gauge-invariant variable 
(density perturbation in synchronous gauge) shown in 
FigJ3] support this conclusion. (We do not discuss the os- 
cilations at early stage visible in Fig. 3 of [12J because of 
different initial conditions and background.) The ratios 
of quintessence (w — const, c\ ~ 1) density perturbations 
to matter density ones in synchronous gauge are shown 
in Fig.l of [|J for fc = O.Olh -1 Mpc -1 . Presented here 
in Fig|5] analogical ratios in conformal-Newtonian gauge 
are similar. The conclusion about anti-correlation be- 
tween the perturbations of the matter and quintessence 
has been done by [17[ and [30( on the base of analy- 
sis of their evolution in the matter rest frame. Q and 
(45| noted this effect too. Recalculation of the frame- 
dependent variables to gauge-invariant ones will - in our 
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FIG. 5: Dependences of ratios Dg'/Dg' on scale factor for 
linear perturbations with scales k — 0.0001, 0.001, 0.01 and 
0.1 Mpc~ x (from top to bottom) in the models with non- 
relativistic matter and quintessence (w = const - top panel, 
c\ — - bottom). 



belief - remove such variance. 



V. CONCLUSION 

The dynamics of expansion of the Universe and 
evolution of scalar perturbations are studied for the 
quintessential scalar fields Q with the classical La- 
grangian L = \Q-iQ' 1 — U{Q) satisfying the additional 
condition w — const or c\ — 0. For both quintessential 
scalar fields the potential U (Q) and time dependence of 
Q are constructed for the same cosmological model and it 
is shown that the accelerated expansion of the Universe 
is caused by the effect of rolling down of the potential to 
minimum (Fig(2J). In QCDM-model with w — const the 
ratio pm/pq —> oo when a — > 0, while in QCDM-model 
with c\ = Pm/pq — > (1 - ^q)/(1 + w )Qq and w — > 
when a — ► 0. At the early epoch w-quintessence is dy- 
namically unsubstantial like cosmological constant while 
c^-quintessence mimics dust matter (w « 0) at a -C 1 
and cosmological constant (w = — 1) at a 3> 1. The de- 
pendence of acceleration parameter on redshift is a bit 
different for them (FigQJ but close to ACDM-model one 
and indistinquishable observationally now. 

Asymptotic analysis of the systems of evolutionary 
equations for gauge-invariant perturbations has shown 
that adiabatic initial conditions for non-relativistic mat- 
ter and w- and c^-quintessence are allowed. The numer- 
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FIG. 6: Evolution of ratios D {M) a init / D^a and */* init 
for linear perturbations with scales k = 0.0001, 0.001, 0.01 
and 0.1 Afpc -1 (from top to bottom) in the models with non- 
relativistic matter and quintessence (c„ = - solid line, w = 
const - dashed line). 



ical integration of these systems give time dependences 
of gauge-invariant variables for matter and quintessence 
scalar perturbations (Figj3l |4]) . The main conclusion de- 
duced from them is following: the magnitudes of the adi- 
abatic matter density perturbations grow like in ACDM- 
model, while for quintessence Dg®\ are constant 
and D^Q\ V^> grow before the particle horizon entry 
but all variables decay after that in such way, that at the 
current epoch they are approximately two orders lower 
than the corresponding quantities for dust matter on su- 
percluster scales. Therefore, on subhorizon scales the 
quintessential scalar field is smoothed out while the mat- 
ter is clustered. 

The quintessential scalar fields studied here suppress 
the growth of matter density perturbations and the 
magnitude of gravitational potential (FigO. In these 
QCDM-models - unlike ACDM ones - such suppression 
is scale dependent and more visible for c^-quintessence. 
Such features of quintessence are important for calcula- 
tions of the matter density power spectrum at different 
redshifts and the power spectrum of CMB temperature 
fluctuations in the range of scales of the late integrated 
Sachs- Wolfe effect. That can be used for interpretation 
of data of current and planned experiments in order to 
identificate the nature of dark energy. 
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